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INTRODUCTION 


Application  of  the  Parabolized  Navi er-Stokes  (PNS)  technique  to  predict 
the  three-dinensi onal  viscous  flow  about  slender  bodies  of  revolution  at 
supersonic  speeds  and  snail  angles  of  attack  has  been  performed  in  the  past 
with  a  considerable  degree  of  success.1  Of  particular  significance  to  the 
U.S.  Army  is  the  capability  that  this  technique  provides  in  the  prediction  of 
the  Magnus  force  and  nonent  for  spinning  shell  at  snail  angles  of  attack.2  3 
Current  efforts  in  the  prediction  of  the  aerodynamics  of  U.S.  Army  projectiles 
in  supersonic  flight  are  directed  towards  extending  the  appl icabi 1 ity  of  the 
PNS  technique  for  the  calculation  of  more  complicated  flows.  Two  such  appli¬ 
cations  are  discussed  in  this  report;  the  first  involves  solution  of  the  flow 
about  a  spinning  axisymetric  projectile  at  six  degrees  angle  of  attack  with  a 
region  of  lee  side  crossflow  separation,  and,  the  second,  the  calculation  of 
the  flow  around  long  L/D  finned  bodies. 

Previous  studies1  2  4  have  presented  accurate  predictions  of  pitch  plane 
aerodynamic  coefficients  for  axisyrnmetric  projectiles  to  angles  of  attack 
greater  than  ten  degrees.  However,  accurate  predictions  of  Magnus  (side) 
force  and  nonent  have  been  limited  to  angles  of  attack  between  zero  and  six 
degrees.  This  is  due,  in  part,  to  the  inability  of  the  Bal dwi n-Lonax  turbu¬ 
lence  nodel 5  to  properly  model  regions  of  lee  side  crossflow  separation. 
Began’  and  Schiff,6  who  applied  the  PNS  technique  to  examine  the  high  angle  of 


' .  W.E.  Sture'-',  D.C.  My  l  in  and  C.C.  Bus':,  "Computational  Parametria  Study  of 
the  Aerodynamics  of  Spinning  Bodies  at  Supersonic  Speeds,"  U.S.  Army 
Ballistic  Research  Laboratory ,  Aberdeen  Proving  Ground,  Maryland,  ARBRL- 
TR-02366 ,  August  1961.  (AD  A106974) 

l.  W.3.  Sturek  and  L.B.  Schiff,  "Computations  of  the  Magnus  Fffect  for 
Slender  Bodies  in  Supersonic  r lou : U.S.  Army  Ballistic  Research 
Laboratory ,  Aberdeen  Proving  Ground,  Maryland,  ARBRL-TR-02384 ,  December 
1961.  (AD*  Al 100 161 

3.  L.B.  Schiff  and  W.B.  Sturek,  "Numerical  Simulation  of  Steady  Supersonic 
Flou'  Over  an  Ogive  Cylinder  Boattail  Body,"  U.S.  Army  Ballistic  Research 
Laboratory,  Aberdeen  Proving  Ground,  Maryland,  ARB.RL-TR-02363  ,  September 
1981.  (AD  A106060 / 


4.  M.  McWerter,  R.W.  Noack  and  W.L.  Oberkampf,  "Evaluation  of  Inviscid' 
Boundary  Layer  and  Parabolized  N avier-Stokes  Solutions  for  Design,  of 
Reentry  Vehicles ,"  AIAA  Paper  No.  84-9486 ,  22nd  Aerospace  Sciences 
Meeting,  January  1984. 


B.S.  Belavin  and  H.  Lomax, 
for  Separated  Turbulent  FI 
Sciences  Meeting,  January 


"Thin  Layer  Approximation  and  Algebraic 
oxs,"  AT  A.-.  Paver  No.  76-267,  16  th  Aerosra 
1976. 


F.  Began:,  and  L.B.  Schiff,  "Computation  of  Supersonic  'Viscous  Flovs  Arcun  i 
Pointed  Bodies  at  Large  Incidence AIAA  Paper  Kc.  83-0034,  21st  Aerospace 
Sciences  Meeting,  January  1983. 


attack  flow  about  cones  and  a  Secant-Ogive-Cylinder  (SOC)  projectile  configur¬ 
ation,  concluded  that  near  regions  of  crossflow  separation  the  values  of  the 
eddy  viscosity  obtained  from  the  Baldwin  and  Lonax  model  were  nuch  too  high. 
To  overcome  turn's  problem,  they  applied,  with  a  significant  degree  o*  success, 
an  approach  which  limits  the  eddy  viscosity  to  values  that  are  more  realis¬ 
tic.  In  this  report,  modifications  similar  to  those  of  Degani  and  Schiff  are 
applied  to  computations  for  a  Secant-Ogi ve-Cyl i nder-Boattai 1  (SOCBT)  projec¬ 
tile  configuration  at  G .3  degrees  angle  of  attack,  Mach  3,  with  and  without 
spin.  The  results  are  compared  to  predictions  obtained  with  the  original 
Baldwin  and  Lomax  turbulence  model  and  to  wind  tunnel  measurements. 

The  second  set  of  results  presented  in  this  report  concerns  the  applica¬ 
tion  of  the  PNS  technique  to  predict  the  aerodynamics  of  long  L/D  finned 
bodies.  Quite  recently,  several  papers  have  reported  the  application  of  the 
PNS  technique  to  examine  the  aerodynamic  of  finned  body  configurations.7’10 
The  computations  for  the  finned  body  configuration  examined  in  Reference  10, 
resembling  a  long  L/D  kinetic  energy  (KE)  penetrator  projectile,  were,  in  the 
opinion  ot  the  authors,  not  fully  successful.  Solutions  over  the  finned  por¬ 
tion  of  the  body  were  only  obtained  using  large  amounts  of  smoothing  and  a 
coarse  grid  near  the  fin  surface.  Additionally,  it  was  later  discovered  that 
increased  circumferential  grid  resolution  is  required  over  the  long  L /P| 
axi  symmetric  forebody.  The  results  presented  in  this  report  address  these 
areas  of  deficiency  in  the  previous  calculation.  Predicted  pitch-plane  aero¬ 
dynamic  coefficients  for  long  L/D  (20-35  caliber)  axisymmetric  kinetic  energy 
penetrator  forebodies  at  Mach  4  and  5,  and  two  degrees  angle  of  attack,  are 
presented  and  compared  with  available  wind  tunnel  data.  Significantly  im¬ 
proved  results  are  then  presented  for  a  complete  finned  body  configuration  at 
Mach  numbers  3,  4  and  5,  and  two  degrees  angle  of  attack. 


II.  COMPUTATIONAL  TECHNIQUE 

Calculation  or  the  flow  field  over  the  body  is  accomplished  using  the 
Parabolized  Navi er-Stok es  technique.  The  PNS  technique  allows  the  solution  to 


7 .  M.V.  Rai  and  D.S.  Chaus  see ,  "New  Implicit  Boundary  Procedures:  Theory 
and  Applications ,"  AIAA  Paper  No.  83-0123,  21st  Aerospace  Sciences 
Meet  in  a,  January  1983  . 

8.  y.M.  Rai,  P.S.  Chaussee  and  Y.M.  Risk,  "Calculation  of  Viscous  Supersonic 
Flows  Over  Finned  Bodies,"  AIAA  Paper  No.  83-1667 ,  Danvers,  MA ,  < July 
1983. 

9.  J.’’.  Jettmar,  "Calculations  of  Viscous  Supersonic  Flow  Over  Finned  Bodies 
"si':  :  a  '  Thin-Fin’  Approximation,"  AIAA  Paper  No.  84-2114,  AIAA 
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he  spacially  narched  along  the  body  in  the  main  flow  direction  due  to  the  par¬ 
abolic  nature  of  the  governing  equations.  An  initial  plane  of  data  is  re¬ 
quired  to  begin  the  space  marching  procedure  and  nay  be  obtained  either  from 
an  auxiliary  calculation  or  from  a  conical  starting  procedure,  as  has  beer 
done  for  the  results  presented  here.  Both  the  space  marching  and  conical 
starting  procedures  are  outlined  below. 

A.  Space  Marching  Procedure 

The  thin-layer  Parabolized  Navier-Stokes  computati onal  technique  devel¬ 
oped  by  Schiff  and  Steger*1  has  been  employed  to  ca^ulate  the  flow  downstream 
of  the  nose  tip.  The  governing  steady  thin-layer  equations  in  strong  conser¬ 
vative  form  and  generalized  coordinates  are  written  below: 

d^s  5F  as  i  as  . 

3s  +  3n  +  3s  *e  3;  ’ 

where  £,  n,  s  are  the  generalized  coordinate  variables  as  displayed  in  Figure 
1  and  defined  below. 


£  =  s(x)  is  the  longitudinal  (marching)  coordinate 
n  =  n(x,y,z)  is  the  circumferential  coordinate 
S  =  ;(x,y,z)  is  the  near  normal  coordinate. 


This  vector  equation  represents  the  thin-layer  approximation  to  the  equa¬ 
tions  of  mass,  momentum  and  energy  conservation  in  the  three  coordinate  direc¬ 
tions.  The  inviscid  flux  vectors  Fs,  F  and  G  and  the  matrix  of  viscous  terms, 

S,  are  functions  of  the  dependent  variables  represented  by  the  vector, 
q(p,  pu,  pv,  pw,  e) ,  where  p  is  the  density;  u,  v  and  w  are  the  velocity  com¬ 
ponents  in  the  three  spacial  directions  x,  y  and  z;  and  e  is  the  total  energy 
per  unit  volume. 

The  parabolized  Navier-Stokes  equations  are  solved  using  a  conservative, 
approximately  factored,  implicit,  finite-difference  numerical  algorithm  as 
formulated  by  Beam  and  Warming.12  Further  details  of  the  numerical  method  may 
be  found  in  Reference  11.  Fitting  of  the  outer  bow  shock  has  been  performed 


21.  L.B.  Scki^f  and  tT.L.  Steaer,  "Nunerical  Si’vulatio n  or  Steads  Supersonic 
Viscous  Flo w,"  A1AA  Paper  No.  79-1131,  17 th  Aerospace  Sciences  'teetinz, 
January  1?79. 

1 F.  R.  Bean-  and  r,.~.  War”:  in  ; ,  "*n  Implicit  Factored  Scheme  r or  the  Compress- 
ibie  Navier-Stokes  Equations AIAA  Journal,  Vol.  IF,  No.  4,  IP7 6  ,  pp. 


in  those  cal cul at  ions,  and  details  of  the  implicit  boundary  procedure  as 
implemented  by  Rai  and  Chaussee  may  be  found  in  Reference  7. 

F.  Con-icai  Star+inn  Solutions 


The  initial  plane  of  data  required  to  begin  the  marching  procedure  is  ob¬ 
tained  US'* ng  the  marching  code  by  assuming  conical  flow  at  the  tip  of  the  pro¬ 
jectile.  By  selecting  a  conical  grid  and  initially  setting  the  flow  field 
variables  to  the  free  stream  values,  the  solution  is  marched  one  step  down  the 
body.  The  solution  is  then  scaled  back  to  the  original  station  according  to 
the  conical  f 1  ow  assumption  and  again  marched  a  single  step.  This  procedure 
is  repeated  until  a  converged  solution  is  obtained.  The  convergence  criterion 
for  the  conical  starting  solutions  applied  here  was  that  the  change  in  density 
between  successive  iterations  was  less  than  10’5  times  the  free  stream  value 
for  each  of  the  points  on  the  body.  This  converged  solution  is  then  used  as 
the  initial  plane  of  data  in  the  marching  procedure. 

It  should  be  noted  that  for  calculations  involving  spinning  projectiles, 
the  conical  starting  procedure  introduces  a  small  error  since  the  circumferen¬ 
tial  velocity  at  the  body  surface  changes  with  longitudinal  position,  viola¬ 
ting  the  conical  flow  assumption.  This  error  is  small,  however,  and  the  cor¬ 
rect  circumferent i a  1  velocity  at  the  body  surface  is  accounted  for  as  the  so¬ 
lution  is  marched  downstream. 

C.  Turbulence  Mode1 


A  fully  turbulent  boundary  layer  has  been  simulated  in  each  of  the  re¬ 
ported  PUS  calculations  using  a  two-layer  eddy  viscosity  model.  For  the  long 
l/D  forebody  and  the  finned  body  configurations,  the  original  turbulence  model 
o*  Baldwin  and  Lomax5  was  applied.  Moderate  angle  of  attack  solutions  for  the 
axi symmetric  shell  configuration  were  obtained  using  both  the  original  Baldwin 
and  Lomax  model  and  the  modified  model.  Both  of  these  models  are  described 
bel ow. 

1.  Baldwin  an-!  Lomax  Model .  The  two-layer,  Cebeci-type,  algebraic  tur¬ 
bulence  model  reported  by  Baldwin  and  Lomax5  been  used  with  considerable  suc¬ 
cess  throughout  many  applications  of  the  thin-layer  Navi er-Stokes  codes.  The 
model  accounts  for  the  effects  of  turbulence  through  an  eddy  viscosity  coeffi¬ 
cient,  and  utilizes  a  purely  local  analysis  of  the  turbulent  region. 

In  the  original  formulation  of  the  Bal dwi n-Lomax  turbulence  model, 
the  turbulent  eddy  viscosity  coefficient,  u^,  is  defined  separately  in  an 

inner  region,  close  to  the  body,  and  an  outer  region,  as  shown  below; 


inner  y  <  y^ 

^ Ut ^  outer  y  <  y 


(2) 


where  y  is  the  radial  distance  from  the  body  surface  and  yc  is  the  smallest 
value  y  at  wh'ck  the  inner  and  outer  values  of  are  equal. 


1? 


Within  the  inner  region,  the  turbulent  eddy  viscosity  coefficient  is 
defined  using  a  nixing  length  node! 


where  the  length  scale,  £,  i 


Whinner 
defined  as 


p 


£2  |U| 


£  =  ky  [1  -  exp(-y+/A+)]. 


(3) 


(4) 


k  and  A+  are  constants,  |w|  is  the  magnitude  of  the  local  vorticity  vector, 
and  y*  is  the  nondinensional  boundary  layer  coordinate  defined  below, 


y 


w  w 


y- 


w 


(5) 


The  subscript  "w"  represents  the  values  at  the  body  surface. 

In  the  outer  region,  for  attached  boundary  layers,  the  turbulent  eddy 
viscosity  is  defined 


outer  =  K  CcppFWAKE  FKLEB(y) 

where  K  and  Ccr)  are  constants,  and  FFLEB(y)  is  the  Klebanoff  internittency 
factor,  which  is  equal  to  1  at  the  wall,  and  decreases  out  from  the  wall. 

^ WAKE  is  3iven  by 


^WAKE  =  ymax  ^"max 


(7) 


where  Fnax  and  ynax  are  determined  from  the  function 


F (y)  =  y  M  [1.  -  exp  (-y+/A+)] 


(8) 


F 


na  x 


is  the  maximum  value  of  F(y)  along  a  radial  profile  from  the  wall  to  the 


boundary  layer  edge.  Tbe  value  of  y  at  which  this  maximum  occurs  is  defined 
as  ynax.  Further  details  of  the  model  can  he  found  in  Reference  5. 

For  analysis  of  the  turbulent  boundary  layer  along  projectiles  at 
small  angles  of  attack  (0  <  a  <  4),  the  Bal dwi n-Lomax  model  performs  satisfac¬ 


torily  due,  in  pari,  to  ine  we 1 1 -oenavea  variations  ot  i-max  and  ymax,  both 

axially  and  circumferentially.  The  lack  of  crossflow  separation  allows  a 
local  profile  of  F(y)  to  attain  a  single,  well-defined  peak.  Accuracy  in  the 
selection  of  the  exact  location  of  a  peak  in  F(y)  is  enhanced  by  fitting  a 
curve  between  three  radial  grid  points  at  every  body  location. 

2.  Turbulence  Model  Modification.  In  order  to  obtain  successful  pre¬ 
dictions  of  flow  fields  containing  regions  of  crossflow  separation,  modifica¬ 
tions  to  the  original  Baldwin  and  Lomax  model  have  been  necessary.  The  ap¬ 
proach  reported  by  Degani  and  Schiff6  has  been  adopted  here  in  an  initial 
attempt  to  compute  the  complex  flow  behavior  associated  with  crossflow  separa¬ 
tion  and  to  determine  the  extent  to  which  such  a  strategy  may  be  applied  to 
spinning  and  nonspinning  boattailed  projectiles. 

In  Reference  6,  computations  for  a  nonspinning  Secant-Ogive-Cylinder 
(no  boattail)  projectile  configuration  revealed  that  the  values  of  the  outer 
eddy  viscosity  calculated  by  the  original  Bal dwi n-Lomax  model  were  much  too 
high  in  the  regions  of  crossflow  separation.  The  cause  was  determined  to  be 
an  ambiguity  in  determining  the  peak  in  F(y)  (Equation  8).  Over  much  of  the 
leeward  side  of  the  projectile,  the  radial  profile  of  F(y)  developed  a  second, 
extraneous  peak  due  to  the  shedding  of  a  vortex  sheet  associated  with  cross- 
flow  separation  at  six  degrees  angle  of  attack  and  greater.  The  behavior  of 
the  function  F(y)  caused  a  meaningless  length  scale  to  be  chosen.  This  sane 
behavior  of  the  function  F(y)  was  also  repeatedly  observed  during  the  course 
of  the  current  study. 

The  strategy  implemented  here  is  similar  to  that  of  Degani  and 

Schiff:  limit  the  values  of  Fnax  and  ymax  obtained  by  the  original  model  by 

anticipating  the  uncontrolled  growth  of  the  function  F(y).  For  each  axial 

station,  a  maximum  value  of  the  scaling  length  ymgx  is  defined  as  1.8  times 
the  value  of  ynax  on  the  windward  ray. 

The  implementation  proceeds  for  each  roll  angle,  starting  from  the 
windward  ray  and  ending  along  the  leeward  ray.  When  the  local  ymax  exceeds 

the  maxinun  allowed  ymax,  it  is  "clipped"  to  the  maximum  allowed  value.  Fmax 

and  ymax  are  then  frozen  with  respect  to  roll  angle  until  the  local  value  of 

ymax  is  again  found  to  lie  within  the  maximum  allowed  value.  A  peak  in  F(y) 

is  defined  if  the  value  of  F(y)  drops  below  90%  of  the  local  maximum.  Where 
two  separate,  distinct  peaks  in  F(y)  exist,  the  peak  closer  to  the  body  is 
chosen.  If  the  two  peaks  in  F(y)  merge  into  one  abnormally  large  peak  (or  a 

peak  cannot  be  found  at  all),  the  value  of  Fmax  is  frozen  at  the  value  used 

fo"-  the  previous  roll  angle  and  the  value  of  ynax  is  set  equal  to  l.s  times 

the  value  on  the  windward  side. 


part, 


the  well-behaved  variations 


Figures  2  and  3  show  the  representat i ve  behavior  of  Fnax  and  ynax 

for  the  original  and  modified  Ba 1 dwi n-Lonax  formulations  for  the  nonspinning 
S0CR7  at  six  degrees  angle  of  attack.  At  each  longitudinal  position,  constant 
values  of  Fmax  and  ymax  are  assigned  over  much  of  the  leeward  face  of  the 

projectile.  The  longitudinal  variation,  although  not  shown  here,  varies 
smoothly  since  it  is  controlled  by  the  smooth  behavior  on  the  windward  side. 
It  is  worth  noting  here  that  the  cutoff  points  for  Fmax  and  ymax  in  Figure  2 

do  not  necessarily  represent  the  roll  angles  where  crossflow  separation  actu¬ 
ally  occurs;  they  merely  reflect  the  factor  of  1.8  being  introduced  into  the 
selection  process.  Furthermore,  interpolating  between  three  radial  points  to 
increase  the  accuracy  of  the  peak  locations  of  F(y)  was  discontinued  where- 
ever  the  cutoff  strategy  was  applied. 


III.  RESULTS 

A.  Shell  Configuration  at  Moderate  Angle  of  Attack. 

Results  were  obtained  primarily  for  the  six  caliber  SUCBT  configuration 
shown  in  Figure  4  at  Mach  3,  C.3  degrees  angle  of  attack,  with  and  without 
spin,  for  fl ow  conditions  duplicating  that  of  the  experiments.13"15  Addition¬ 
al  results  we>-e  also  computed  for  the  same  projectile  with  no  boattail  (SOC) 
to  establish  effects  of  presence  of  a  7°  boattail. 

The  grid  for  these  computations  consisted  of  45  exponentially  stretched 
points  between  the  body  and  the  shock  and  an  equal  spacing  of  points  circum¬ 
ferentially  at  5°  increments.  For  the  case  of  zero  spin,  calculations  were 
performed  in  a  half  plane  by  applying  a  symmetry  boundary  condition  along  the 
midplane  of  the  body,  while  the  computations  with  spin  were  carried  out  in  a 
full  circumferential  plane.  To  properly  resolve  the  viscous  effects,  the  ra¬ 
dial  grid  resolution  at  the  body  surface  was  controlled  such  that  the  nondi- 
mensional  boundary  layer  coordinate,  y+,  was  approximately  5  at  the  first 
point  above  the  wall  for  all  roll  angles.  This  was  accomplished  through  the 


13.  R.P.  Reklis  and  W.B.  Sturek ,  "Surface  Pressure  Measurements  on  Slender 
Bodies  at  Angle  of  Attack  at  Supersonic  Speeds,"  U.S.  Army  Ballistic 
Research  Laboratory,  Aberdeen  Proving  Ground,  Maryland,  ARBRL-MR-02B76 , 
November  1978.  (AD  A064097 ) 

14.  C.J.  Nietubicz  and  K.O.  Opalka,  "Supersonic  Wind  Tunnel  Measurements  of 
Static  and  Magnus  Aerodynamic  Coefficients  for  Projectile  Shapes  with 
Tangent  and  Secant  Ogive  Noses,"  U.S.  Army  Ballistic  Research  Laboratory , 
Aberdeen  Provina  Ground,  Maryland,  ARBRL-MR-02991 ,  February  1980.  (AD 

AO 8329? ) 

15.  L.D.  Kayser  and  W.B.  Sturek,  "Turbulent  Boundary  Layer  Measurements  cm. 
the  Boattail  Section  of  a  Yawed,  Spinning  Projectile  Shape  at  ’-tick  3.0," 
U.S.  Army  Ballistic  Research  Laboratory,  Aberdeen  Proving  Ground, 
Maryland,  ARBRL-MR-02880 ,  November  1978.  (AD  A06S355) 
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adaptive  grid  technique  first  reported  in  Reference  2.  The  computational  grid 
was  generated  al gebraical ly  for  this  axisynnetric  configuration. 


The  marching  step  size  and  smoothing  parameters  employed  for  the  SOC  and 
SOCBT  configurations  were  chosen  according  to  the  results  presented  in  Refer¬ 
ence  10.  The  marching  step  size  at  the  nose  tip  was  .01  calibers  and  was 
increased  by  5%  every  15  steps  to  a  value  of  .03  calibers  at  the  tail  end  of 
the  projectile.  Only  fourth-order  explicit  smoothing  was  applied  in  these 
computations  and  was  set  equal  to  5  times  the  step  size.  All  SOC  and  SOCBT 
results  were  generated  on  a  CDC  7600  computer  with  a  speed  of  .00247  CPU 
sec/step/grid  point.  The  total  run  times  were  approximately  30  minutes  for  a 
nonspin  case  and  60  minutes  for  a  spin  case. 

Computations  were  performed  by  applying  either  the  original  Bal dwi n-Lomax 
turbulence  model  or  the  modified  turbulence  model  over  the  entire  length  of 
the  projectile. 

1.  Nonspin  Case.  For  the  case  of  zero  spin,  the  initial  appearance  and 
downstream  growth  of  the  separated  region  on  the  lee  side  of  the  SOC  and  SOCBT 
projectiles  can  be  clearly  traced  in  Figures  5a  and  5b.  The  circumferential 
location  of  crossflow  separation  is  defined  as  the  roll  angle(s)  at  which  the 
circumferential  shear  at  the  wall  changes  sign.  The  region  of  crossflow  sepa¬ 
ration  appears  on  the  leeward  side  at  an  axial  location  of  approximately  4 
calibers  and  immediately  begins  to  widen  circumferentially  around  the  projec¬ 
tile.  The  predicted  axial  position  where  crossflow  separation  initially 
appears  and  the  size  of  the  separated  region  is  somewhat  sensitive  to  whether 
the  original  or  modified  model  is  applied.  It  is  also  apparent  that  a  major 
effect  of  the  boattail  is  to  cause  this  separated  region  to  broaden  further 
around  the  projectile. 

The  computed  longitudinal  surface  pressure  distributions  on  the  wind 
and  lee  sides  of  the  SOCBt  at  a  =  6.3°  are  compared  to  experiment  in  Figure 
6.  The  modified  turbulence  model  was  used  for  these  computations.  The  agree¬ 
ment  is  quite  good,  most  notaoly  in  the  vicinities  of  discontinuities  in 
stream  wise  surface  curvature.  The  well-behaved  longitudinal  variation  of 
surface  pressures  is  used  here  as  an  indicator  of  the  suitability  of  the 
marching  step  size  and  the  smoothing  parameters  chosen. 

The  computed  circumferential  surface  pressure  distributions  upstream 
of  the  boattail  agree  well  with  experiment  with  or  without  the  modifications 
to  the  turbulence  model,  as  was  previously  reported  in  Reference  6.  In  con¬ 
trast,  Figure  7  shows  a  comparison  of  circumferential  surface  pressure  distri¬ 
bution  near  the  end  of  the  boattail  configuration  for  the  original  and  modi¬ 
fied  turbulence  models  compared  with  experimental  measurement.  The  improve¬ 
ment  in  agreement  due  to  the  modified  turbulence  model  is  twofold.  First,  the 
agreement  in  trend  on  the  leeward  side  is  much  improved;  and  second,  the  mag¬ 
nitude  of  the  pressure  on  the  windward  side  agrees  within  1%. 

The  general  effect  of  the  modified  turbulence  model  on  longitudinal 
velocity  profiles  is  shown  in  Figure  8.  For  zero  spin,  at  roll  angles  where 
the  turbulence  scaling  lengths  are  "clipped",  a  distinct  deficit  in  the  pro¬ 
file  is  consistently  generated. 


For  zero  spin,  Figures  9a,  9b  and  9c  show  longitudinal  velocity  pro¬ 
files  near  the  end  of  the  boattai  1  for  the  original  and  modified  turbulence 
models  compared  to  experiment.  These  three  roll  angles  are  included  in  the 
region  of  separation  as  defined  in  Figure  5b.  At  roll  angles  120°  and  150°, 
the  modified  turbulence  model  gives  acceptable  agreement  with  experiment  and 
shows  favorable  trends  compared  to  the  original  turbulence  model.  At  roll 
angle  180°,  however,  the  modified  model  gives  a  less  favorable  agreement  with 
experiment.  This  implies  that  the  modified  turbulence  model,  as  implemented 
here,  does  not  provide  an  appropriate  length  scale  near  roll  angle  180°. 

Very  little  difference  in  the  predicted  aerodynamic  forces  and 
moments  was  observed  using  either  the  original  or  modified  turbulence  model 
for  the  nonspinning  shell. 

2.  Spin  Case.  A  more  sensitive  test  of  the  performance  of  the  modified 
turbulence  model  is  the  prediction  of  the  Magnus  effect  for  a  spinning  shell. 
Application  of  the  modified  turbulence  model  to  the  spin  case,  it  was  hoped, 
would  follow  directly  from  the  case  of  zero  spin.  However,  initial  results 
indicated  that  the  cutoff  distance  for  ymax  (1.8  times  the  value  on  the  wind¬ 
ward  side)  was  too  strict  for  the  case  of  spin.  By  relaxing  the  cutoff  dis¬ 
tance  to  3.0  times  the  value  on  the  windward  side,  better  agreement  was  ob¬ 
tained  between  computation  and  experiment. 

The  development  of  Magnus  force  (sign  convention  given  in  Figure  10) 
over  the  spinning  SOCBT  projectile  is  shown  in  Figure  11.  Both  the  original 
and  nodified  turbulence  models  yield  good  agreement  with  experiment.  If,  when 
applying  the  modified  model,  the  cutoff  distance  factor  for  ynax  had  been  set 

to  1.8  rather  than  3.0,  the  Magnus  force  would  have  been  overpredicted  by  20%. 

The  effect  of  the  modified  turbulence  model  on  the  longitudinal  ve¬ 
locity  profiles  for  the  spin  case  compared  to  experiment  is  shown  in  Figures 
12a  thru  12c.  As  for  the  nonspin  case,  the  nodified  turbulence  model  causes  a 
velocity  deficit  compared  to  the  original  model.  At  roll  angle  210°,  the 
modified  model  agrees  substantially  better  with  experiment  than  the  original 
model.  At  roll  angles  120°  and  240°,  the  overall  agreement  is  satisfactory. 
At  roll  angle  180°  (not  shov/n  here)  the  agreement  is  less  favorable. 

B.  Long  L/D  Finned  Body  Configurations 

The  results  for  long  L/D  finned  body  conf igurations  are  presented  in  two 
parts.  First,  results  are  presented  for  smooth  axisymmetric  forebodies  of  a 
shape  typical  of  kinetic  energy  penetrators  and  comparisons  made  with  wind 
tunnel  data.  Results  are  then  presented  for  a  complete  finned  body  configura¬ 
tion  (forebody  plus  fins)  and  comparisons  made  with  inviscid  code  predictions 
and  with  free- flight  range  data. 

1.  Forebodies.  Computations  have  been  performed  using  the  PUS  technique 
to  examine  the  effect  of  length  to  diameter  ratio  (L/D)  on  the  aerodynamics  of 
axisymmetric  kinetic  energy  penetrator  forebodies  and  the  results  compared 
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with  wind  tunnel  data. 16  Figure  13  shows  the  forebody  shape  as  modeled  in  the 
wind  tunnel  tests.  The  primary  difference  between  the  computational  model  and 
the  wind  tunnel  model  is  that  the  small  step  decrease  in  diameter  over  the  aft 
portion  of  the  wind  tunnel  model  forebody  has  been  neglected  in  the  computa¬ 
tional  model.  This  difference,  it  is  felt,  should  not  have  a  significant 
effect  on  the  prediction  of  pitch-plane  aerodynamic  coefficients. 

The  qridding  for  these  calculations  consisted  of  60  points  from  the 
body  to  the  shock  and  37  points  circumferenti ally  around  the  body  in  the  half 
plane  of  symmetry.  The  ci rcumferent i al  gridding  was  chosen  after  performing 
a  study  examining  the  effect  of  ci rcumferent ial  resolution  on  an  eight  degree 
cone-cylinder  configuration,  for  free-flight  conditions,  Mach  4,  and  free 
stream  Reynold's  number  (based  on  body  diameter)  of  3.2  million.  Figure  14 
shows  the  development  of  the  normal  force  over  the  body  for  varying  degrees  of 
ci rcumferenti a1  resol ution.  While  19  circumferential  points  may  be  adequate 
for  calculations  of  short  L/D  bodies,  such  as  shell  configurations,  long  L/D 
bodies  require  more  resolution  for  accurate  solution.  Circumferential  resolu¬ 
tion  of  317  points  and  above  appears  to  give  adequate  resolution  for  these  long 
axisymmetric  bodies. 

Results  are  presented  here  for  Mach  numbers  of  4  and  5,  2C  angle  of 
attac1',  and  turbulent  flow  conditions  over  the  body.  Free  stream  conditions 
correspondi ng  to  the  wind  tunnel  test  conditions  have  been  used.  Constant 
wall  temperature  boundary  conditions,  (wall  temperature  =  294K)  were  also 
applied  at  the  body  surface. 

After  obtaining  the  initial  plane  of  data  from  the  conical  marching 
procedure  at  a  position  0.53  calibers  from  the  nose,  the  solution  was  marched 
down  the  body  using  a  step  size  of  .021  calibers,  The  spacing  from  the  wall 
to  the  first  point  above  the  wall  was  adjusted  so  that  the  first  point  above 
the  wall  was  within  the  laminar  sublayer;  the  boundary  layer  coordinate,  v+ , 
was  maintained  less  than  five. 

Figure  15  and  16  show  the  development  of  the  normal  force  and  pitch¬ 
ing  moment  (referenced  to  the  projectile  nose)  over  the  body  at  Mach  4  and  2° 
angle  cf  attack.  Good  agreement  is  seen  between  the  computation  and  the  ex¬ 
periment  at  25  calibers.  Scatter  in  the  experimental  data  at  35  calibers  is 
evident,  emphasizing  the  difficulties  in  performing  such  a  test  for  these  long 
bodi es. 

Similar  agreement  is  seen  for  the  results  at  Mach  5  and  2°  angle  of 
attack,  shown  in  Figures  17  and  18.  Again,  scatter  in  the  experimental  data 
is  evident  at  the  higher  L/D 1  s . 

Figure  19  shows  the  effect  of  Mach  number  and  Reynold's  number  on 
the  development  of  normal  force  over  the  body.  Here,  the  results  discussed 
above  for  Mach  4  (Re  =  0.41  million,  based  on  body  diameter)  and  Mach  5  (Re  = 
0.54  million)  are  shown  along  with  a  free-flight  calculation  at  Mach  4  (Re  = 


1C.  F.  Brandon,  "private  communications,"  unpublished  wind  tunnel  data,  l 
Arty  Ballistic  Research  Laboratory ,  Aberdeen  Brovina  Ground ,  Marxian . 
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2.2  million).  Evidently  ->  this  Mach  number  and  Reynold's  number  regime,  both 
Mach  number  and  Reynold's  number  are  significant  parameters  to  be  simulated. 
While  computations  can,  co  some  degree,  simulate  both  Reynold's  and  Mach  num¬ 
bers,  such  simulation  is  difficult  in  the  wind  tunnel  and  suggests  the  need  to 
scale  the  wind  tunnel  data  to  free  flight  conditions  for  these  long  L/D 
bodi es. 

The  effect  of  length  to  diameter  ratio  on  the  development  of  normal 
force  can  also  be  deduced  from  Figure  19.  For  shorter  bodies,  L/D  <  12,  the 
development  of  normal  force  increases  significantly  with  L/D,  and  the  vari¬ 
ation  is  due  almost  purely  to  inviscid  effects.  For  larger  L/D's,  the  vari¬ 
ation  of  normal  force  with  L/D  is  smaller,  and  is  dominated  primarily  by 
viscous  effects.  For  length  to  diameter  ratios  greater  than  about  20,  the 
normal  force,  in  the  range  of  Reynold's  numbers  of  interest,  shows  a  small 
increase  due  to  the  thickening  of  the  boundary  layer. 

2.  Finned  Projectile  Configuration.  The  finned  body  configuration  for 
which  calculations  have  been  performed  closely  resembles  the  M735  Army  projec¬ 
tile.  The  modeled  finned  body  configuration  is  characterized  by  a  conical 
nose  section  joined  tc  a  smooth  ey7  •  nd^i cal  main  body  with  six  symmetrical 
swept  fins  attached  to  the  aft  section  of  the  projectile.  Figure  20  depicts 
the  basic  dimensions  of  this  configuration.  The  actual  projectile  differs 
from  the  modeled  projectile  in  that  the  actual  projectile  has:  (1)  circumfer¬ 
ential  grooves  over  much  of  the  cylindrical  portion  of  the  body  to  prevent  the 
sabot  from  sliding  off  the  body  in  the  gun  tube;  (2)  fins  which  have  a  non- 
synmetrical  sectional  geometry  to  induce  roll,  and  (3)  a  slightly  rounded 
nose.  Modeling  the  projectile  with  a  sharp  nose  and  symmetrical  fin  section 
are  not  the  result  of  inherent  limitations  of  the  computational  model,  but 
rather  a  matter  of  convenience  for  these  initial  calculations.  While  modeling 
of  the  sabot  grooves  may  be  possible  using  surface  blowing,  wind  tunnel  re¬ 
sults  have  shown  that  such  grooves  have  almost  no  effect  on  the  value  of  nor¬ 
mal  force  and  pitching  moment.  These  grooves  do,  however,  have  a  noticeable 
effect  on  drag,  particularly  at  higher  angle  of  attack.  16 

A  shadowgraph  of  the  actual  projectile  in  flight  at  Mach  4.3  is 
shown  in  Figure  21  and  displays  some  of  the  relevant  features  of  the  flow 

field;  a  bow  shock  wave  emanating  from  the  nose  of  the  projectile,  shocks  at 
the  leading  edge  of  the  fins,  expansion  waves  at  the  cone-cylinder  junction, 
and  a  boundary  layer  which  increases  along  the  body. 

Results  are  presented  here  for  Mach  number  of  3,  4  and  5,  2°  angle 
of  attack,  and  turbulent  flow  conditions  over  the  body.  Atmospheric  flight 
conditions  were  simulated  by  maintaining  the  body  temperature  at  the  free 

stream  value  of  294 K.  Calculations  were  made  with  two  of  the  fins  oriented 
vertically,  enabling  a  half  plane  of  symmetry  to  be  applied. 

An  initial  plane  of  data  was  generated  using  the  conical  step  back 
procedure  at  position  0.35  calibers  from  the  nose.  The  solution  was  then 

marched  down  the  body  to  a  position  near  the  beginning  of  the  fins  using  a 

step  size  of  0.014  calibers,  60  points  from  the  body  to  the  shock,  and  121 
points  circumferential ly  about  the  body  in  the  half  plane  of  symmetry.  Over 
the  cone-cylinder  portion  of  the  body  only  fourth-order  explicit  snooth'nq  was 
used.  Spacing  from  the  wall  to  the  first  point  above  the  wall  was  again  main¬ 
tained  so  that  at  least  one  point  was  in  the  laminar  sublayer. 


Once  on  the  finned  portion  of  the  body,  the  grid  was  obtained 
through  the  use  or  an  elliptic  grid  generator.7  Points  on  the  body  surface 
were  clustered  near  the  leading  edge  of  the  fins,  as  shown  in  Figure  22. 
Figures  23a  and  b  show  a  cross  section  of  the  grid  on  the  finned  portion  of 
the  body  at  an  axial  location  on  X/D  =  13.2.  Grid  points  are  clustered  near 
the  body  to  resolve  the  boundary  layer. 

Substantial  difficulty  had  been  encountered  in  previous  calcula¬ 
tions*0  in  marching  the  solution  over  the  finned  portion  of  the  body.  Prior 
solutions  were  only  obtained  using  large  amounts  of  smoothing  arid  a  coarse 
grid  which  did  not  properly  resolve  the  viscous  layer  over  the  fins.  For  the 
current  set  of  calculations,  a  small  fillet  was  applied  axially  fron  the  cone- 
cylinder  onto  the  leading  edge  of  the  fins.  The  fillet,  of  constant  radius, 
began  0.07  calibers  in  front  of  the  fins  and  ended  at  a  position  where  it  be¬ 
came  tangent  to  the  leading  edge  of  the  fins.  This  allowed  a  solution  to  be 
obtained  with  significantly  improved  resolution  of  the  boundary  layer  over  the 
fins  and  with  only  5%  of  the  previously  applied  smoothing.  Current  values  of 
smoothing17  used  for  these  calculations  are  shown  in  Table  1.  In  marching  the 
solution  over  the  finned  portion  of  the  body,  the  step  size  was  reduced  from 
0.014  calibers  to  0.0072  calibers. 


Table  1.  Smoothing  Parameters  for  Finned  Body  Calculations 


SMOOTHING 

PARAMETERS  CONE-CYLINDER 


FI  MS 

M  =  3  M  =  4  M  =  5 


SMI’ 

0.015 

0.05 

0.05 

0.025 

SM'JIM 

0. 

0.1 

0.1 

0.05 

EPSA 

0. 

0.05 

0.05 

0.025 

CP  SB 

0. 

0.05 

0.05 

0.025 

Figure  24  displays  the  pressure  along  the  body  on  the  wind  and  lee 
sides  and  at  a  roll  angle  of  90°  at  Mach  4.  A  sharp  rise  in  pressure  on  the 
leading  edge  of  the  wind  and  lee  fins  is  seen,  followed  by  a  sharp  drop  at  the 
axial  location  where  the  fin  reaches  its  maximum  span.  Over  the  finned  por¬ 
tion  of  the  body  the  4  =  90°  line  shows  a  jump  in  pressure  due  to  the  inter¬ 
action  of  the  shocks  from  the  adjacent  fins.  It  is  interesting  to  note  that, 
due  to  the  significant  reduction  in  smoothing  required  to  generate  a  solution, 
the  pressure  predicted  on  the  leading  edge  of  the  fins  is  smaller  than  that 


Numerical 
it  m-/>:  A. 
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previously  calculated,  and  the  jump  in  pressure  due  to  the  interaction  of  the 
shocks  is  more  pronounced. 


Figure  25  shows  the  development  of  the  normal  force  coefficient  over 
the  body.  The  normal  force  coefficient  shows  a  moderate  contribution  due  to 
the  conical  nose  and  cylindrical  portion  of  the  body  and  a  larger  contribution 
due  to  the  finned  portion  of  the  body.  In  this  figure  comparison  is  made  with 
results  obtained  with  the  NSWC  inviscid  code  (SWINT)18  and  range  data19  for 
the  actual  projectile  configuration.  The  SWINT  code  solves  the  Euler  equa¬ 
tions  for  the  supersonic  flow  over  bodies  with  fins  and/or  wings  and  for  the 
external  flow  about  bodies  with  inlets.  The  code  makes  use  of  the  thin  fin 
approximation,  collapsing  each  fin  along  a  single  radial  plane  in  the  grid. 
The  fin  leading  edges  must  be  sharp  and  cannot  extend  beyond  the  bow  shock. 
Variations  in  the  fin  cross  section  geometry  can  be  accounted  for  in  the  code 
by  the  application  of  appropriate  local  analysis  such  as  shock  compression  and 
Prandtl-Meyer  expansion  theories,  though  this  has  not  been  utilized  in  the 
computations  presented  here.  For  the  current  SWINT  computations,  the  fins 
have  been  modeled  as  a  flat  plate  of  zero  thickness. 

Development  of  the  normal  force  as  predicted  by  the  PNS  and  SUIN'!' 
codes  compare  well,  particularly  over  the  cone-cylinder  portion  of  the  body. 
For  forebodies  with  L/D's  greater  than  the  current  configuration,  differences 
in  the  viscous  (PNS)  and  inviscid  (SWINT)  code  predictions  can  be  expected  due 
to  viscous  effects.  (Figure  19,  previously  presented,  shows  little  effect  of 
Reynold's  number  at  L/D's  less  than  about  12.)  Over  the  finned  portion  of  the 
body,  slightly  more  lift  is  being  predicted  in  PNS  computations.  Good  agree¬ 
ment  between  the  total  value  of  normal  force  for  both  procedures  is  seen  com¬ 
pared  with  the  range  data. 

The  variation  of  the  predicted  values  of  the  slope  of  the  normal 
force  and  pitching  moment  coefficients  with  Mach  number  are  compared  with  the 
values  predicted  by  the  SWINT  code  and  with  range  data  in  Figures  26  and  27. 
(Pitching  moment  here  has  been  referenced  to  the  center  of  gravity  shown  in 
Figure  20.)  The  PNS  code  is  seen  to  predict  slightly  larger  values  of  normal 
force  and  pitching  moment  compared  with  the  range  data  and  inviscid  code  pre¬ 
dictions  at  Mach  3  and  4.  PNS  and  inviscid  code  predictions  agree  very  well 
at  Mach  5.  Differences  in  agreement  between  the  PNS  and  inviscid  code  predic¬ 
tions  at  Mach  3  and  4  occur  almost  entirely  over  the  finned  portion  of  the 
body,  though  at  this  point  such  differences  can  not  be  attributed  entirely  to 
viscous  effects.  It  should  be  noted  again  that  the  exact  fin  geometry  has  not 
been  modeled  in  either  the  PNS  or  SWINT  computations.  Efforts  to  model  the 
fin  geometry  as  closely  as  possible  are  underway. 
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An  additional  calculation  at  Mach  3  was  run  over  the  finned  portion 
of  the  body  using  131  points  circumferentially  in  the  half  plane.  Variations 
in  the  predicted  values  of  normal  force  between  the  computations  with  121  and 
181  ci rcumferential  points  was  less  than  a  tenth  of  a  percent.  Comparisons  of 
the  circumferential  pressure  distribution  at  several  axial  stations  showed 
very  little  difference.  This  result  established  that  the  current  circumferen¬ 
tial  resolution  of  121  points  was  adequate. 

IV.  CONCLUSIONS 

A.  Computations  at  Moderate  Angle  of  Attack 

Computational  results  for  SOC  and  SOCBT  shell  configurations  with  and 
without  spin  have  been  generated  at  angle  of  attack  6.3°.  One  set  of  results 
employs  a  modification  to  the  Baldwin-Lomax  turbulence  model;  the  other  set  of 
results  uses  the  model  in  its  original  form.  Both  sets  of  results  were  com¬ 
pared  to  available  wind  tunnel  data  in  order  to  assess  the  predictive  capa¬ 
bility  for  flows  with  crossflow  separation  at  this  angle  of  attack. 

For  the  case  of  zero  spin,  the  modified  turbulence  model  significantly 
improved  the  agreement  of  ci rcumferenti al  pressure  distributions  where  cross- 
flow  separation  is  present.  Velocity  profile  comparisons  on  the  boat tail  were 
generally  improved,  with  the  modified  model  consistently  generating  a  velocity 
deficit  compared  to  the  original  model.  The  variation  between  the  predicted 
pitch  plane  and  drag  coefficients  using  the  original  and  modified  turbulence 
models  was  very  small,  and  the  coefficients  were  in  good  agreement  with 
experiment. 

For  the  case  with  spin  (PD /V  =  .19),  the  modified  turbulence  model  had  to 
be  relaxed  to  obtain  good  agreement  with  experimental  Magnus  force  and  veloc¬ 
ity  profile  data.  The  need  to  allow  larger  scaling  lengths  for  the  case  of 
spin  is  evidence  that  characteristic  differences  may  exist  in  the  turbulence 
modeling  requirements  for  spinning  and  nonspinning  bodies  with  crossflow 
separati on. 

B.  Computations  for  Long  L/D  Finned  Body  Configurations 

Results  of  the  computations  for  the  long  L/D  axi symmetric  forebodies  show 
good  agreement  with  the  wind  tunnel  data  for  pitch-plane  aerodynamic  coeffi¬ 
cients  at  L/D's  of  20  and  25.  Scatter  in  the  experimental  data  at  higher 
L/D's  makes  comparisons  difficult.  The  computations  demonstrate  the  impor¬ 
tance  of  considering  viscous  effects  for  these  long  bodies. 

Computations  were  performed  for  the  finned  body  configuration  and  show 
encouraging  agreement  with  range  data  and  inviscid  code  predictions.  Computa¬ 
tions  were  obtained  with  significantly  improved  resolution  of  the  viscous 
layer  on  the  fins  and  without  the  need  for  large  amounts  of  smoothing  over 
this  portion  of  the  body,  compared  to  results  previously  obtained. 
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Figure  2.  Circumferential  Variation  of  Fmax  for  S0C  and  S0CBT 

Configurations;  M  =  3,  a  =  6.3°,  PD /V  =  0, 
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Figure  Longitudinal  Surface  Pressure  Distribution  for  SOCBT 
Configuration;  M  =  3,  a  =  6.3°,  PD/V  =  0, 
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Figure  7.  Circumferential  Surface  Pressure  Distribution  for  SOCBT 
Configuration;  M  =  3,  a  =  6.3°,  PD/V  =  0, 

Re  =  2.13  x  107/m,  X/D  =  5.77 
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Figure  3.  Longitudinal  Velocity  Profiles  for  SOCBT  Configuration;  M 
a  =  6.3°,  PD/V  =  0,  Re  =  2.13  x  107/m,  X/D  =  5.67 
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Figure  9a.  Longitudinal  Velocity  Profile  for  SOCBT  Configuration;  M 
a  =  6.3°,  PD/V  =  0,  Re  =  2.13  x  107/m,  X/D  =  5.67,  a  = 
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Figure  12a.  Longitudinal  Velocity  Profile  for  SOCBT  Configuration;  K  =  3, 
a  =  6.3°,  PD/V  =  .19,  Re  =  2.13  x  107/n,  X/D  =  5.67,  *  =  120° 
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Figure  12c.  Longitudinal  Velocity  Profile  for  SOCBT  Configuration;  M  = 
a  =  6.3°,  PD /V  -  .19,  Re  =  2.13  x  107/m,  X/D  =  5.67,  9  =  2 
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Figure  13.  Forebody  Configuration 
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Figure  19.  Development  of  Normal  Force  Coefficient  over  Long  L/D  Axi symmetric 

Body,  Effect  of  Mach  and  Reynold's  Numbers,  a  =  2° 
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speed  of  sound 

pitching  moment  coefficient 

dCm/da,  slope  of  the  pitching  moment  coefficient 
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dC^/da,  slope  of  the  normal  force  coefficient 

Magnus  (side)  force  coefficient 
diameter  of  model 

total  energy  per  unit  volume  of  fluid,  normalized  by  p^a^ 

flux  vectors  of  transformed  gas  dynamic  equation 
projectile  length 
Mach  number 

pressure  normalized  by  p^a^ 

nondimensicnal  spin  rate  about  model  axis 
Reynold's  number  per  unit  length,  P^M^a^/p^ 

Reynold's  number,  p  a  D/p 

00  00  CO 

Reynold's  number  based  on  diameter  p  M  a  D/u 
viscous  flux  vector 

Cartesian  velocity  components  along  the  x,  y,  z  axis, 
respectively,  normalized  by  a^ 

physical  Cartesian  coordinates 

distance  from  wall  in  law  of  the  wall  coordinates 

angle  of  attack 

coefficient  of  viscosity,  normalized  by  free  stream  value 
turbulent  eddy  viscosity  coefficient 
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to 


density,  normalized  by  free-stream  density 
shear  stress 

circumferential  angular  coordinate 
local  vorticity  vector 
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•  =  free-stream  conditions 

w  =  body  surface  values 
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